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, Abstract 

(N 

' The influence of nuclear multiscattering and shadowing on pion spectra is investigated in pA 

I collisions from CERN SPS to RHIC energies. The calculations are performed in a next-to-leading 

order (NLO) pQCD-improved parton model, including intrinsic partonic transverse momentum 
I distributions. The nuclear modification of the pion spectra (Cronin effect) is considered at 

I different targets in a wide energy range. Theoretical predictions are displayed for planned pA 

experiments at CERN SPS and recent dAu experimental results are analysed at RHIC. 



O ; 1 Introduction 

The Cronin effect was discovered in the energy range 20 GeV < ^/s < 40 GeV at FERMILAB P 
121 O]- The effect is a ~ 20 — 40% enhancement in particle production in proton-nucleus (pA) 
^ ' collisions compared to scaled proton-proton (pp) ones, where the scaling factor is the number of 

binary collisions. Recent RHIC experiments at ^/s = 200 GeV raised the interest on this nuclear 
effect in nucleus-nucleus (AA) collisions. Theoretical investigations were launched to understand 
and reproduce the Cronin peak [SElini- 

Theoretical models were based on perturbative Quantum Cromodynamics (pQCD). General 
problem with pQCD models, that they underestimate the experimental data by ~ 50 — 100% at 
usual scales. Changing scales could solve this problem, however larger scale causes an increase in the 
lower limit of applicability of pQCD calculations over the region of interest 3 GeV< < 6 GeV. 
Introducing the so called intrinsic transverse momenta (/ct)i one can reproduce the experimentally 
measured pion (vr) spectra in pp collisions within ~ 10 — 20%. The relevant parameter of these 
calculations is (k^), the width of the intrinsic partonic transverse momentum distribution which 
depends on the cm. energy {^/s), at fix scales. 

Here pion production is investigated in a wide (high-)energy range from CERN SPS up to RHIC 
(20 GeV < ^/s < 200 GeV). Although this analysis can be revisited in case of other hadrons|ni Ej, 
however due to the poor statistic of experimental data, in this paper we would like to reproduce 
only pionic data. Our calculation is performed in next-to- leading order (NLO), and compared to 
earlier leading order (LO) results 0. 



2 Theoretical Background for pp collisions 



First of all let us start with the invariant cross section for pion production in a pp collision, which 
can be described in the fc^-enhanced NLO pQCD-improved parton model on the basis of the 
factorization theorem as a convolution [H]: 

""d^p^ S^^Jvw/z^ v{l-v) Jvw/vzc w J 



^6{l-w) + Kab,c{s, V, W, Q, Qr, Qf) 

dv n 



(1) 



where we introduced factorised 3-dimensional parton distribution functions (PDFs), 

/(x,kT,Q') = f{x,Q^)-g{kT) . (2) 

Here, the function f{x,Q'^) represents the standard 1-dimensional NLO PDF as a function of mo- 
mentum fraction of the incoming parton x at factorization scale Q, da/dv represents the Born cross 
section of the partonic subprocess ah — > cd, Kab,c{s,v,w,Q,QR,QF) is the corresponding higher 
order correction term, and the fragmentation function (FF), D'^{zc,Q'f), gives the probability for 
parton c to fragment into a pion with momentum fraction Zc at fragmentation scale Qf- We use 
the conventional proton level {S, V, W) and parton level (s, v, w) kinematical variables of NLO 
calculations (see Ref.s 01101 )■ 

In this analysis we consider fixed scales: the factorization and the renormalization scales are 
connected to the momentum of the intermediate jet, Q = Qr = k ■ pq (where pq = pt/zc), while 
the fragmentation scale is connected to the final hadron momentum, Qf = k ■ Pt- The value of k 
is ~ 0(1). 

For our NLO calculations we are using a 3-dimensional PDF which is a product of the standard 
PDF (at its Q scale) and a 2-dimensional initial transverse-momentum distribution, g{kT) of par- 
tons containing its "intrinsic kx" parameter as in Refs. [HllHlHiEj- We demonstrated the success 
of such a treatment at LO level in Ref. 6 , and a Kjet-hased NLO calculations in jlll I12j . In our 
phenomenological approach the transverse-momentum distribution is described by a Gaussian, 

ff(kT) = ^e-'^T/C^^) . (3) 

Here, {k'^) is the 2-dimensional width of the kx distribution and it is related to the magnitude of 
the average transverse momentum of a parton as (kj^) = 4(fcT)^/7r. 

In present article we use the LO GRVd^, and the NLO MRST(cg)[Tl] PDFs. For FF we apply 
most recent KKP parameterizations ^3], both in LO and NLO cases. These PDF and FF sets can 
be applied down to very small scales (Q^ « 1.25 GeV^), thus we have got results at relatively small 
transverse momenta, ^ 2 GeV at our fix scales. The details of the NLO calculations can be 
found in Ref. 

In our earlier work|S] we determined the pT-independent intristic-Zcr values for the 3 GeV 
< < 6 GeV momentum window and investigated their dependence on cm. energy at fixed 
scales [TT]. We found that the reproduction of the Cronin effect requires (kj^) ~ 2 GeV^. The di-jet 
production data at ISR energies yield a similar value for {k'^) ^0]- Recent measurements of jet-jet 
correlations in pp collision at RHIC energies |17l I18j can clarify the properties of the transverse 
component of the PDFs. 

Our description of pion production in pp collision can be extended for pA and AA collisions. 



3 Nuclear Effects in pA and AA Collisions 



Proton-nucleus and nucleus-nucleus collisions can be described by including collision geometry, 
saturation in nucleon-nucleon {NN) collision number, and shadowing inside the nucleus. In the 
frame of Glauber picture, the cross section of pion production in nucleus-nucleus collision can be 
written as an integral over impact parameter b: 



where pp cross section on the right hand side represents the cross section from eq. but with 
an increased widths compared to the original transverse-momentum distributions ^ in pp colli- 
sions, as a consequence of nuclear multiscattering (see eq. ©). Here tj^{b) = J dzpj^{b,z) is the 
nuclear thickness function (in terms of the density distribution of nucleus A, pa), normalized as 
/ d'^btAib) = A. For small size nucleus we used sharp sphere approximation, while for larger nuclei 
Wood-Saxon formula were applied. 

The initial state broadening of the incoming parton's distribution function is accounted for by 
an increase in the width of the Gaussian parton transverse momentum distribution in eq. ©: 

(^t)pA = {kT)pp + C ■ hpA{b) . (5) 

Here, {kj)pp is the width of the transverse momentum distribution of partons in pp collisions, hpAib) 
describes the number of effective NN collisions at impact parameter b, which impart an average 
transverse momentum squared C. The effectivity function hpA{b) can be written in terms of the 
number of collisions suffered by the incoming proton in the target nucleus, i^Aib) = (7NNtA{b), 
where is the inelastic NN cross section: 

hpAib) = \ ^^^'^ . (6) 
I J^m — 1 otherwise 

The value Vm = oc corresponds to the case where all possible semihard collisions contribute to the 
broadening. We have found that for realistic nuclei the maximum number of semihard collisions is 
3 < ^^„^ < 4. 

The determination of the factor C and of VA{b) is in progress in a systematic analysis in 
NLO jni- Our preliminary results confirm the findings of Ref. jSj; where the systematic analysis 
of pA reactions was performed in LO and the characteristics of the Cronin effect were determined 
at LO level. Following Ref. IS], we assume that only a limited number of semi-hard collisions (with 

= 4) contributes to the broadening, average momentum squared imparted per collision is C = 
0.4 GeV2. 

Furthermore, the PDFs are modified in the nuclear environment by the "shadowing" effect [2111 
1211 122j . This effect and isospin asymmetry are taken into account on average using a scale inde- 
pendent parameterization of the shadowing function SajA{x) adopted from Ref. [1]: 



fa/A{x,Q^) = Sa/A{x) 



jfa/pix,Q^)+(l-j)fa/nix,Q'') 



(7) 



where fa/n{x,Q'^) is the PDF for the neutron and Z is the number of protons. In the present 
work, we display results obtained with the EKS parameterization [2J, and with the updated HI- 
JING parameterization |2()j . The former has an antishadowing feature, while the later incorporates 
different quark and gluon shadowing, and has an impact-parameter dependent and an impact- 
parameter independent version. The impact-parameter dependence is taken into account by a term 
oc (1 — b'^/R\), which re- weighs the shadowing effect inside the nucleus. 



We note, that for the deuteron (d), one could use a simple superposition of a pAu and a nAu 
collision, or a distribution for the nucleons inside the deuteron. For a first orientation we follow 
Ref. [231 124j in this regard, and apply a hard-sphere approximation for the deuteron with A = 2 
taking no nuclear effects for the deuteron. 

4 Results on the Cronin Effect 

Including these effects one can calculate the invariant cross section for an AA' collision. Thus, 
introducing the nuclear modification factor Raa' ) as 

_ da^^' I d^p ( "included nuclear effects" ) 
dfj^^'/d^P ("i^o nuclear effect") ' 

nuclear effects can be investigated clearly and efficiently, on a linear scale. 

Following^ 12] for the nuclear modification factor, in Figure 1, we present the Cronin effect as 
a ratio of pW to pBe collisions {RwiBe)- This "type" of nuclear modification factor, is principally 
insensitive to the shadowing at CERN SPS energies. We compared experimental data^ 12113], with 
theoretical calculations at ^fs = 19.4, 23.8, 27.4 and 38.8 AGeV cm. energies. In all calculation we 
choose {kj') ~ 2 GeV^, 1/^ = 4 and Csat = 0.4 GeV^ for all A, while the scale was fixed as k = 1/2 
in case of LO calculations {dotted lines). In NLO case the scales are chosen as «; = 2/3. It seems 
there is no difference (within ~ 5%) between curves with no shadowing {solid lines) and with the 
updated HIJING 6-independent {dashed) or 6-dependent {dash-dotted) lines. 

Figure 2 presents the dependence of RpA {left panel) and Ra/Bc {right panel) on different targets 
A G { ^Be, 32^, 44^^, 98/n and ^^'^Au } at CERN SPS energies in NLO calculations at Ebeam = 158 
AGeV, with (kj,) = 1.9 GeV^ . We present our results with EKS shadowing {dotted lines), the 
6- independent, and the 6-dependent updated HIJING shadowings {solid lines). The 5-independent 
and 5-dependent updated HIJING cases are coincide within errors. 

In our model the position of Cronin peak does not depend explicitely on the cm. energy, 
however, depends on the value of (kj-). With the above scale choices the peak is appears at ~ 4 
GeV. Furthermore, the height of the peak depends on the A via C ■ hpA{b) term in eq. El (see 
Fig. 1). In the energy range 20 GeV < ^/s < 40 GeV shadowing effects are generally small. Thus, 
the influence of the different shadowings on the Cronin effect are negligible, especially in minimum 
bias cases (see Fig. 2). 

In Figure 3 we display our NLO results for the nuclear modification factor in dAu collision 
for 6-independent {solid line) and 6-dependent {dashed line) shadowing parameterizations from the 
updated HIJING shadowing using parameters {^/s = 200 AGeV, (kj.) = 2.5 GeV^, k = 4/3). In 
spite of the very different impact parameter dependence of RdAu the minimum bias results are 
very close to each other. This is because many details are averaged out in minimum bias data, 
and the final result is no longer sensitive to the 6-dependence of shadowing. As we can see on 
Fig. 3 the theoretical results reproduce quite well the minimum bias experimental data on vr^ from 
PHENIX|25j. In minimum bias case we do not see a strong Cronin peak. However, a reasonable 
Cronin peak will appear in different centrality bins. 

Figure 4 summarizes our results on the impact-parameter (6) dependence of the Cronin effect 
in dAu collision |24j. The nuclear modification factor, RdAu{b), is presented with updated HIJING 
shadowing. The left column shows the 6- independent cases, while the right column displays the 
6-dependent ones. 

Applying 6-independent shadowing for central collisions {left upper panel) the curves overlap, 
since in these bins, the target matter is approximately constant from the point of view of multiscat- 
tering and shadowing. In more peripheral cases {left lower panel), shadowing starts to dominate, 
RdAu{b) is decreasing from bin to bin. The 6-independent shadowing suppresses RdAu{b) under 1. 



(8) 



In case of 6-dependent shadowing at central bins {right upper panel) Cronin effect wins against 
suppression. At peripheral bins (right lower panel) - where Cronin effect begin to vanish and 
shadowing suppression working also -, the nuclear modification factor becomes unity. 

5 Summary 

We investigated pp and pA experimental data on pion production with our /cT-augmcntcd pQCD 
parton model at NLO level in the cm. energy window 20 AGeV < ^/s < 200 AGeV. In the lower 
limit of this energy range (-y/s < 40 AGeV), the pion p^-spectra is not sensitive to the applied 
shadowing mechanism, because shadowing is close to be negligible (~ 5%) in this kinematical 
region. However the multiscattering strongly influences this spectra and dominates every deviation 
of pA collision from superimposed nucleon-nucleon collisions We determined the most important 
parameters of the parton model from existing experimental data and we predicted the Cronin effect 
in pA collisions at CERN SPS energy {Ehgam = 158 AGeV). Considering different targets {^Be, 
32g^ 44^1^^ ^^/n and ^^'''Au) we obtained an increasing multiscattering effect with increasing size of 
the target nucleus, while the position of the Cronin peak did not change. 

Our analysis was repeated at RHIC energy {y/s = 200 AGeV). Here the shadowing effect is much 
larger, thus the differences between the competing shadowing models are emphasized. Applying 
6-dcpcndent shadowing description for the dAu collision the centrality dependence of the Cronin 
effect generates an unusual structure in the nuclear modiflcation factor, which can be measured 
experimentally. In this way we can distinguish between different shadowing models from the data. 
However minimum bias calculations yields similar results for different shadowing models. Using 
our model we successfully reproduced the measured nuclear modification factor for 7r° in the dAu 
collision. 
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FIG. 1. Cronin effect on Rw/Be nuclear modification factor at CERN SPS energies, from ^/s — 
19.3 AGeV up to 38.8 AGeV compared to experimental dataQlEllSl- Curves are in LO with dotted, 
in NLO with no shadowing (SO) with solid, and in NLO with updated HIJING shadowingj^OI both 
in 6-independent dashed and 6-dependent dash-dotted hnes cases. The later two are not different 
within errors. All calculations were carried out using = 4, Csat = 0.4 GeV^, and with {k'^) 
shown. We applied scales: k = 1/2 in LO and k = 2/3 in NLO. 
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FIG. 2. Dependence of Raa^Vt) {left panel) and RA/BeiPr) {right panel) on different A targets 
{A e {^Be, 325, 44Ti, 98 j„ and ^^T^u}) at CERN SPS energies. We presented our results on 
three types of sliadowing: EKS by Eskola et aLj^lj {dotted lines), and new HIJING!^ shadowing 
{solid lines). The NLO calculations were carried out at Eteam = 158 AGeV, (k^) = 1.9 GeV^ and 
at K = 2/3. 
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FIG. 3. The nuclear modification factor RdAu for minimum bias dAu collisions obtained from our 
NLO calculation in case of updated HIJING 6-independent, {solid line) and 6-dependent shadowing 
{dashed line). Theoretical calculations are compared to two different experimantal sets PbGI {dots) 
and PbSc {open circles) at PHENIX 
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FIG. 4. The impact parameter (b) dependence of the nuclear modification factor RdAu{b) with 
5-independent shadowing (right column) and with 6-dependent shadowing [left column) . The fines 
correspond to different centrafity bins. The 3 "central" curves (0-50%) are shown in the upper 
panels, and the remaining 3 "peripheral" cases (50-100%) can be seen in the lower panels. 



